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— Random Lattice Fermions at E = without Doubling — 
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Random bond Hamiltonians of the it flux state on the square lattice are investigated. It has 
a special symmetry and all states are paired except the ones with zero energy. Because of this, 
there are always zero-energy states (or zero-modes). The states near the zero-modes are described 
by massless Dirac fermions. For the zero-mode in a system with odd numbers of sites, we can 
construct a random lattice fermion without a doubling and quite large systems ( up to 801 x 801) 
are treated numerically. We clearly demonstrate that the zero-mode is given by a critical wave 
function. For the multifractal behavior of the wave function, the agreement with the effective field 
theory is demonstrated. 

PACS numbers: 73.40H, 61-43-j, 72-15R 



Random critical points in two dimensions appear in 
quantum Hall systems and in systems with spin-orbit 
coupling. Both critical points have important experimen- 
tal implications, and both are in strong coupling limit. 
Because of that we are still unable to calculate their crit- 
ical exponents analytically, despite the efforts of more 
than a decade. Recently, a class of exactly solvable two- 
dimensional random critical points was found. |,§ Those 
critical points appear in models of Dirac fermions coupled 
to random gauge fields. However, it was pointed out that 
those critical points contain infinite number of relevant 
directions. Thus it seems that it is almost impossible to 
realize those random critical points in a lattice model. 

In this paper, we study random bond models on the 
square lattice with only nearest neighbor hopping. The 
models have n flux per plaquette on the average. When 
all hopping matrix elements are real (in this case the 
model will be called random hopping (RH) model), we 
find that with weak randomness and in a continuum ap- 
proximation, the states near zero energy are described by 
two copies of Dirac fermions coupled to imaginary ran- 
dom chemical potential, real random mass, and imagi- 
nary gauge fields. The results in Ref. suggest that 
the continuum model has a line of random critical points, 
and due to the imaginary gauge fields and the imagi- 
nary chemical potential, the critical points has only one 
marginally relevant direction if we fix the energy to be 
zero. Thus it could be possible to fine tune parameters 
to reach the critical line. In fact, the zero-energy states 
of the fine tuned model are found to have a multifractal 
behavior and are critical. Q] We also present direct nu- 
merical results indicating that the zero-energy states in 
the generic random bond model also appear to be criti- 
cal with multifractal behaviors at least for the systems of 
sizes up to 801 x 801. The multifractal scaling function 



f(a) for the fine-tuned RH model is calculated numeri- 
cally which agrees with the exact result of the continuum 
theory. GIB . Another interesting point in our treatment 
is that we can avoid the fermion doubling as far as the 
zero-energy states are concerned. 

In one dimension, some of the wave functions of the 
quasiperiodic systems show the critical behavior clearly. 
H In two dimensions, it was difficult to show the critical- 
ity clearly due to the limitation of available system sizes. 
However, since our models have several specialties as de- 
scribed below, we can treat quite large two dimensional 
systems up to 801 x 801 which enables a clear demon- 
stration. 

The wave functions in a two-dimensional random sys- 
tem are believed to be always localized if the system has 
the time reversal symmetry. Thus it is quite interest- 
ing that our models have random critical points. There 
are several reasons why the critical states are allowed in 
our models in two dimensions. Our Random bond model 
has a very special property that the random Hamiltonian 
anti-commute with the operator 7 defined below. This 
means that the eigenstates always appear in pairs with 
energies E and —E. Therefore the states near zero en- 
ergy are quite special, and one may wonder if they are 
critical states or not. We may also consider a pure sys- 
tem with a small diagonal hopping m (which break the 
symmetry mentioned above) . |^| The symmetry is recov- 
ered by taking a limit m — > 0. When m ^= 0, there is an 
energy gap near E = and there are two energy bands. 
The Hall conductance of the lower band is +sign(m) and 
— sign(m) for the higher one. Since the nonzero Hall con- 
ductance is carried by extended states, there are paired 
extended states in the lower and the higher bands. By 
taking m — > 0, the paired extended states merges and 
disappear at E = 0. It suggests the criticality of the 
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E = states. Also from a point of the perturbation the- 
ory, the density of states of the pure state vanishes at 
zero energy. Therefore a usual treatment of the purta- 
bative consideration has to be modified. This property 
we believe, is also responsible for the appearance of the 
critical states at zero energy. 
Our Hamiltonian is 



H 



E 

<i,j> 



+ h.c. 



(1) 



where the summation is over the nearest-neighbor bonds. 
We shall analyze the following two types of random 
bond models near the it flux: (i) random gauge model 
t j+x<j = (-)3v e M*U) } t j+yd = e m y^\ (ii) random hop- 
ping (RH) model, tj +x j = (— ) ly + St x , tj +y j = 1 + St y . 
Here 59 XtV (j) and 6t x y are real random variables. Nu- 
merical results of the wave functions for these two mod- 
els are qualitatively similar. The model (i) belongs to 
the unitary ensembles. The model (ii) preserves the time 
reversal symmetry and belongs to the orthogonal ensem- 
bles which is suitable for a comparison with the critical 
continuum theory 

In the absence of the randomness, the system is invari- 
ant under the translation and the Hamiltonian is written 
in the momentum space: 



cos ky cos k% \ 
'■ \ cosfc, -cosfc y I ^ k ( > 



with ip\ — (c£, c£ + / s). We see that the zero-energy 
state is formed by the fermions with momentum k\ = 
(7r/2,7r/2) and k 2 = (— tt/2, tt/2). Expanding it near k\ 
and fc 2 , we find in the continuum approximation, H be- 
comes 



H = 2iJ 



d 2 x& 











<J 3 

CT 3 



(3) 



with *t = 1/(4tt 2 ) / d5k x 5k y e- lSk * x ~ lSk y y {ipi i7 ipl 2 ) and 
Cj x _j y — (p x+ iy , 0*-iv , i-Ox+jy ^ i-jx-jy )\|> where cr's are 
the Pauli matrices. We assume that take the same 
value on four sites (x, y) — (2m, 2n), (2m + 1, 2n), 
(2m, 2n + 1), (2m + 1, 2n + 1). Note that for weak im- 
purities, the states near zero energy are described by the 
smooth function $?(x). Therefore we can derive a contin- 
uum theory to describe those states near zero energy. 

Before deriving the continuum theory with random- 
ness, we first study the symmetries of the lattice Hamil- 
tonian. Since the square lattice is bipartite, the trans- 
formation c jxJy -> c' jxjjy = {-l) ix+jy Cj xJv induces 
sign change of the Hamiltonian. It is expressed as 
H(x) — > 7 t iJ(x)7 = —H(x), 7 = oi <g> o\ where 
H = J d 2 x^(x)^ H(x) 1 i(x) is the continuum Hamilto- 
nian. The contribution from the randomness SH = 
H — Ho has a form 5H = J d 2 x^g(x)^ if we ignore the 



terms containing derivatives. The above condition im- 
plies that the most general form of the g(x) is given by 
a linear combination of the 4 by 4 matrices: 71 = o 2 <8> I, 
72 = o~i ® <t 2 j 73 = — <T 2 ® o~i, 74 = I ® o~i with coeffi- 
cients 2cii(x), i = 1, ...,4. The total Hamiltonian, after 
a change of basis to make 7 = — 03 ® I, ( ^ = T*', 
T = (7 <g> I + ia 2 ® cri)/V2), is given by 



H 



j d 2 x^2(I( 



(4) 



L>- 
D - 

D± = i[(d x =F ai)ai =F (d y T 02)0-3] ± m 3 =F o- 2 a 4 . (5) 



Note that the averaged Hamiltonian is invariant under 
translation (x, y) — > (x+1, y). This requires that (at) = 
for z = 1, 2,3. The averaged Hamiltonian is also invari- 
ant under reflection (x,y) — ► (x,—y). This implies that 
(04) = 0. 

The zero-energy state of if satisfies 



L> + V>+ = 0, £>_V- = 0, 



(6) 



where 



'J''. Thus the random bond model model 



contains two Dirac fermions in a continuum approxima- 
tion, which is the famous fermion doubling of the lattice 
fermions. It causes several difficulties in numerical calcu- 
lations. However, as we discuss below, we can avoid the 
doubling as far as the E = state is concerned. 

Now let us discuss about the doubling which arises 
from the lattice models. Since our model is a nearest 
neighbor hopping model on a bipartite lattice, we can 
write the Hamiltonian as 



H=({cl},{cl}) 



O V 
2?t O 



(7) 



where {c + } = {c jx j y \j x + j y is even } = {c ++ , c — } is a 
set of fermion operators at one of the sublattices (+) and 
{c-} = {c?W« \jx+j y is odd } = {c+_, c_+} is the other 
(— ). Here let us assume a system to be L x x L y with both 
odd L x and L y with a fixed boundary condition. ( We 
label the sites as (1, 1), • • • , (L x , L y ).) Then a total num- 
ber of sites L x L y = N + + iV_ is odd with N+ = N- + 1 
where N + = #{+ sites}, iV_ = #{— sites}. Therefore 
I? is a 7V + x iV_ = N + x (N + — 1) rectangular matrix. 
The Schrodinger equation is reduced to TJD^(t) + = E 2 <f> + 
and X>tX>0_ = E 2 4>-. Any state with E ^ ks 
a pair state at — E as *(</> + , </>_)_£ =* (</> + ,—^>_)e- 
O V* 



Since det 



( note that I? is a form of 



N + x N + — 1), there is always zero-mode. The wave 
function at E = satisfies 



V^l =° = 0, Vcj) 



E=0 



0. 



(8) 
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Since dim(/)^ =0 = iV + , the first one gives JV_ 
equations for iV + variables and there is a one dimen- 
sional non vanishing solution. (Therefore this zero-mode 
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is not degenerate in general.) On the other hand, since 
dim i ? =o = iV_ = N + — 1, the second one gives N + equa- 
tions for iV_ = N + — 1 variables, this is overdetermined. 
Therefore we have <^ =0 = 0. 

It gives an important restriction in the continuum 
model. It implies 

tp- = (9) 

in (|^) since the transformation matrix between |\E')"s 
and the site operators' UT U , only mixes operators 
within each sublattices. Therefore there is no degeneracy 
( doubling ) for the E = state. This E — state is de- 
scribed by a usual two component spinor. Therefore the 
independent component is not four but two as far as the 
E = state is concerned. We could avoid the doubling 
by using the above special conditions. 

Random Gauge Model. — Now let us show numeri- 
cal results for the random gauge model. In the previous 
numerical calculations for the 7r flux state with random- 
ness H , it was not possible to treat large two-dimensional 
systems. In the present work, however, we only need the 
eigenstate of a semi positive definite operator T>T>\ with 
zero eigen value. This and the fact that the E = eigen 
states are generally not degenerate enable us to treat 
large size systems up to 801 x 801. We take 56 x (i) and 
S0 y (i) to be uniform random numbers between —V/2 to 
V/2. In Fig. |, the wave function is shown for V = 0.2 
with 801 x 801 system where it is coarse grained over 4 
sites on the plaquette. It seems to be neither localized 
nor extended in a usual manner. To understand the na- 
ture of the wave functions, the function f(a) obtained 
numerically are also shown. The box size dependence is 
carefully included in the calculations. |5| The maximum 
of f(a) gives a Hausdorf dimension which is always two 
here since the wave functions are on the square lattice. 
The value ao which gives the maximum of f(a) gives 
scaling of the dominant parts of the wave function. It 
is two for an extended states which gives uniform non- 
singular scaling. For a critical state, ao ^ 2 and has a 
singular scaling. In our case, it apparently deviates 
from 2 which means the zero-mode is critical. 

Random Hopping Model. — The properties of random 
Dirac fermion models have been studied in detail lately. 
We study the random hopping model to make a compar- 
ison with the continuum theory. In the continuum ap- 
proximation, the fact that 8t Xty 's are real further requires 
19{ x )l = d{ x )* ■ Thus we take with the coefficients ai(x), 
i = 1, ...,4 to be real. Note that D± happens to be the 
Hamiltonian of Dirac fermions with imaginary random 
gauge potentials (a%, 02), an imaginary random chemical 
potential 04, and a real random mass 03. Let (01,02), 
203 and 2a4 to have a Gaussian distribution with width 
gA, gv and Qm respectively. Averaging over randomness 
generates four-fermion interactions with (gA, gv, 9m) as 
coupling constants. The one-loop beta functions for the 
above coupling constants is calculated in Ref. S as 



g A = 32g M gv (10) 
gv + 9m = ~8(gv ~ 9m) 2 (11) 
9v - 9m = -8(gv + 9m)(9v - 9m) ~ 8g A (gv - 9m) (12) 

(Note comparing to Ref. here we have an imagi- 
nary random chemical potential, a real random mass and 
an imaginary random gauge potential. Thus the beta 
functions are modified accordingly with (17^4 , gv , 5m ) — > 
(— gA, — 9v, 9m)-) The renormalization-group (RG) 
equation has a fixed line (gA,gv,9M) = (gA,®, 0). 
Around the fixed line, gv~9M is irrelevent (note gA > 0) 
and thus we may set gv — 9m- In this case gv + 9M does 
not flow, and non-zero gv+9M will drive gA — > 00. Thus 
the fixed line has at least one (marginally) relevent direc- 
tion. Two-fermion operators are also relevent. But those 
two- fermion operators do not appear since < ai > = 0, 
i = 1, 4 due to the discrete symmetries. Because of the 
imaginary random gauge potential, the results in Refs. 

|H imply that all charge neutral operators are ir- 
relevant, except two-fermion operators and four-fermion 
operators discussed above. Thus the fixed line has only 
one marginally relevent direction. It suggests that the 
RH model have a line of critical points at zero energy 
which is described by the line of critical points in the ran- 
dom Dirac fermion model. Those results further suggest 
that, once we fix the energy to be zero, there is only one 
marginally relevant direction for the critical line. Thus 
the critical line in the RH model contains two relevant di- 
rections. Due to its relation to the random Dirac fermion 
model, the critical line in the RH model is exactly solv- 
able, and critical exponents can be calculated exactly. 

To describe the critical line, set 03 = 04 = 0. Let us as- 
sume the probability distribution of (ai(r), 02 (r)) is given 
by P cx exp(— J d 2 r^—[a\(r) + a|(r)]) where gA, char- 
acterizing the strengtn of the gauge field fluctuations, is 
exactly marginal. Thus gA can be used to parameterize 
different critical points on the critical line. In general 
critical exponents are functions of gA ■ 

Note that the energy E couples to an operator \E , t_\E , _ ( _. 
Since \& + and belong to two independent Dirac 
fermion models described by D + and £>_ respectively, 
The scaling dimension of is simply the sum of the 

scaling dimensions of ^+ and ^ + . From Ref. [[l0| , we see 
that the scaling dimensions of \&+ and are equal to 
1/2, which implies that the scaling dimension of is 
1. This means that the dynamics exponent z = 1 and the 
density of states scales as N(E) cx \E\; both exponents 
are independent of gA- 

In the continuum theory (|^), the zero-energy state 
can be solved exactly [Q. Introducing <f> and \ through 
a n — d^cf) + e^uduX- One can show that the zero-energy 
state (with, for example, 7 = 1) is given by 

oc e^ +x (13) 
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and the probability distribution of <f> and \ is given 
by P oc exp(-/d 2 r^ I [(5^(r)) 2 + [(%(r)) 2 ]).The wave 
function (|13|) with this distribution was studied in de- 
tail in RcfTy §, and was shown to have a multifractal 
structure. The exact multifractal scaling function f(a) 
is defined on d_ < a < d + and is given by 

[d + — a)(a — d_) 



/(a) - 8 



(14) 



where d± = 2(1 ± Vlf) for 9a < 2vr and d + = 
8-<J~p-, d- = for g^4 > 27r. We see that for weak 
randomness (gA < 27r) f(a) is peaked at 



r.,9A 

a = 2 H 

7T 



(15) 



The parameters <Zj(r) in the continuum model is di- 
rectly related to the lattice randomness. We find that 



03 (J) 



= -AStjj+^-y-, a 2 (j) = -ASt jJ+1) (-) i ' +jv , 



A6t jt j +S: (-y<° +j v and a 4 (j) 



-ASt 



jj+y 



with a dimensional parameter A. Therefore if we chose 



St 



%] to satisfy 



St 



-St 



+x+v- Then g v 



- Stj + y j-f-y + x and Stj j-^y — 

9m = and the lattice model 



will be on the critical line. We used this parametriza- 
tion and performed a numerical calculation on the above 
fine tuned RH model. We set St x ,y as Gaussian random 
numbers with a variance V 2 {gA — A 2 V 2 ). In Fig. ||, 
we have plotted a singular scaling function f{a). In the 
calculation, gA is obtained from a value of «o- The f(a) 
obtained from numerical calculation agrees well with the 
exact theoretical result Eq. (|l4|). 

In summary, we find that in the continuum approxima- 
tion, the RH model at zero energy has a critical line with 
only one marginally relevent direction. This analytical 
result is supported by numerical calculations on the fine- 
tuned RH model. Therefore the RH model is described 
by the continuum random Dirac fermions as far as critical 
behavior is concerned. We also find that the zero-energy 
states of the general random bond models also appear to 
be critical according to our numerical results. However, 
the numerical calculations can not eliminate a possibility 
of localization with huge correlation length ( more than 
800 lattice spacings) which may appear in our problem 
due to the marginally relevent direction. 
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FIG. 1. The zero- mode wave function of the random 
gauge model and the corresponding f(a). The system size 
is 801 x 801 and V = 0.2. The errors by the finite size cor- 
rections in the box sizes is about the symbol sizes within the 
shown region but is larger in the other. 



FIG. 2. Functions /(q)'s for the zero- mode wave func- 
tions of the RH model: 401 x 401 and V = 0.3. The solid 
line is the analytical result. The different symbols are for the 
different randomness realizations. 
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